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Abstract 

Subradiance, i.e. the cooperative inhibition of spontaneous emission by destructive interatomic 
interference, can be reahzed in a cold atomic sample confined in a ring cavity and lightened by a 
two-frequency laser. The atoms, scattering the photons of the two laser fields into the cavity-mode, 
recoil and change their momentum. Under proper conditions the atomic initial momentum state 
and the first two momentum recoil states form a three-level degenerate cascade. A stationary 
subradiant state is obtained after that the scattered photons have left the cavity, leaving the 
atoms in a coherent superposition of the three collective momentum states. After a semiclassical 
description of the process, we calculate the quantum subradiant state and its Wigner function. 
Anti-bunching and quantum correlations between the three atomic modes of the subradiant state 
are demonstrated. 

PACS numbers: 03.75.-b, 42.50.Nn, 37.10.Vz 
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I. INTRODUCTION 



Recent experiments with Bose-Einstein Condensates (BEC) driven by a far off-resonant 
laser beam liave demonstrated collective Superradiant Rayleigh and Raman 

scattering {5,0], sharing strong analogies with the superradiant emission from excited two- 



level atoms In these experiments an elongated BEC scatters the pump photons into 
the end-fire modes along the major dimensions of the condensate, acquiring a momentum 
multiple of the two-photon recoil momentum hq, where q = k — kg and k and ks are the 
wave vectors of the pump and the scattered field. Theoretical works have shown that the 
Superradiant Rayleigh Scattering relies on the quantum collective atomic recoil (QCARL) 
gain mechanism, in which the fast escape of the emitted radiation from the active medium 
leads to the superradiant emission 0, yjll^. The quantum regime of CARL 11, 12| occurs 



when the two-photon recoil frequency uj^ = hq^/2m is larger than the gain bandwidth, 
such that the recoil frequency shifts the atoms out of resonance inhibiting further scattering 
processes. As a consequence in the QCARL each atom coherently scatters a single pump 
photon, changing momentum by hq. The process in which the atoms make a transition 
between two momentum states = and p = hq) has strong analogies with that of two- 
level atoms prepared in the excited state and decaying to the lower state by spontaneous 
and stimulated emission. However, the incoherent spontaneous emission dominating the 
two-level atomic decay is absent in the momentum transition, where spontaneous emission 
is associated to momentum diffusion due to the scattering force, which can be made very 
small if the laser is sufficiently detuned from the atomic resonance. The absence of Doppler 
broadening and the long decoherence time of a BEC allows to observe superradiance and 
coherent spontaneous emission much more easily than from electronic transitions in excited 
atoms, in which the decay is dominated by the incoherent spontaneous emission. 

Another example of cooperative phenomena from excited two-level atoms is subradiance, 
i.e. the cooperative inhibition of spontaneous emission by a destructive interatomic inter- 
ference. This phenomenon, whose existence has been proposed by Dicke (1954) in the same 
a.t,c.e predicting supe„adia„ce Q, Kas .ece.ved .esj co^.d^ation t.an t.e more popular 
superradiance, also due to the difficulty of its experimental observation. In fact, the only ex- 
perimental evidence has been done on 1985 by Pavolini et al. 1^. Among different schemes 
of multi- level systems in which subradiance was predicted, Crubellier et al., in a series of 
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FIG. 1: (a): Schematic diagram illustrating the geometry of a two- frequency pump CARL experi- 
ment; (b): three-level cascade scheme. 



theoretical papers proposed a three-level degenerate cascade configuration in 



which cooperative spontaneous emission is expected to exhibit new and striking subradiance 
effects. 

In this paper we show that subradiance in a three-level degenerate cascade can be realized 
in a BEC inserted in a ring cavity and lightened by two laser fields with frequency difference 
twice the two-photon recoil frequency, as illustrated by figiU The frequency of the scattered 
photon is determined by energy and momentum conservation. The process consists in two 
steps. In the first step the atoms initially at rest scatter the laser photons of frequency uj 
into the cavity mode of frequency ujg = u — ujr, changing momentum from to p = hq. 
In the second step the atom scatters the laser photon of frequency u + A changing their 
momentum from p = hq to p ' = 2hq. Since the change of the kinetic energy of the atom 
is AE = [p'"^ — p^)/2m = Shur, by energy conservation the frequency of the scattered 
photon is + A — 3ujr which coincides with the frequency generated in the first step when 
A = 2ujr. In this way a three- momentum- level degenerate cascade is realized in which the 
atoms, initially with momentum p = 0, change momentum to the intermediate value hq 
and then to the final value 2hq, emitting two degenerate photons of frequency Ug = u — Ur- 
In general the process, as described in ref. 19|], will continue with an other scattering of 
the photon of frequency u + A changing the atomic momentum from 2hq to 3hq with the 
emission of a photon of frequency u + A — = uj — Su;^ and so on. However, if the cavity 
linewidth is much narrower than the frequency difference, k <^ 2uj^, these other frequencies 
will be damped out. Then, the oscillation of only the frequency ujs = uo — uor in the cavity 
will restrict the momentum cascade to the three momentum states, 0, hq and 2hq. A basic 
feature of this system is that the transition rates are proportional to the pump intensities, 
so that they can be varied with continuity. This makes the subradiance observation much 
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easier than with a three-level cascade between electronic energy levels, where the transition 
rates are fixed by the branching ratios. 



II. SEMICLASSICAL TREATMENT 



A. General model 



The quantum collective atomic recoil laser (QCARL) with a two-frequency pump is de- 
scribed by the following equations for the order parameter of the matter field and 
the cavity mode field amplitude a{t) fiol : 

,(9* h 



dt 2m dz^ 



+ ig [a(t)a*e*(^"-^*) - cc] ^ (1) 

^ = gNait) [ dz\^\^e'^'^'~^'^ - na, (2) 
at J 

where z is the coordinate along the cavity axis and a{t) = l-|-eexp(— zAt). These equations 
have been derived performing the adiabatic elimination of the atomic internal degrees of 
freedom but replacing the pump field with Ep = e*(*^'^~'^*) (^Eq + Eie"*^*). In Eqs.(Il]) and 
([2]) a{t) = {eoV/2JTLUsy^'^Es{t) is the dimensionless electric field amplitude of the scattered 
radiation beam with frequency tOs, g = {^lo/2Ao){uj(f /2heoVy^'^ is the coupling constant, 
Qq = (IEq/H is the Rabi frequency of the pump laser incident with an angle (p with respect 
to the z axis (0 = vr if counterpropagating) , with electric field Eq and frequency u detuned 
from the atomic resonance frequency ujq hj Aq = u — ujq. The pump laser has a sideband 
with frequency u; + A, with A = 2u!r [where Ur = hq^ /2m and q = 2k sin(0/2)], and electric 
field El with e = Ei/Eq. The other parameters are: d = i ■ d, the electric dipole moment 
of the atom along the polarization direction e of the laser, V, the cavity mode volume, 
the total number of atoms in the condensate, 6 = u — Us, and k, the cavity linewidth. The 
emitted frequency is within the cavity frequency linewidth, whereas the pump field is 
external to the cavity so that its frequencies are not dependent on the cavity ones. The 
order parameter \E' of the matter field is normalized such that / = 1. 

If the condensate is much longer than the radiation wavelength and approximately ho- 
mogeneous, then periodic boundary conditions can be applied on the atomic sample and 
the order parameter can be written as = J2n^n(t)un{z)e~''"'^^, where Un{z) = 

(g/27r)^/^ exp[m(g2;)] are the momentum eigenstates with eigenvalues Pz = n{hq). Using 
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this expansion, Eqs.([T]) and ([2]) become: 

dc 



dt 

where ujn = n{nujr — 5). 



+ g[a{t)a*Cm-i - a*{t)acm+i] (3) 
^ = gNa{t)'^CmC*^^i- na, (4) 



B. Three-level approximation 

As as been discussed elsewhere js], , if the gain rate is smaller than the recoil frequency 
the atoms recoil only along the positive direction of q, absorbing a photon from the laser 
and emitting it into the cavity mode. Backward recoil, in which an atom absorbs a photon 
from the cavity mode and emits it into the laser mode, is inhibited by energy conservation. 
In this way, the laser photon of frequency u induces a momentum transition from m = 
to m = 1, emitting in the cavity a photon with frequency Ug = u — Ur] the laser photon of 
frequency u + A [with A = 2uJr] induces a momentum transition from m = 1 and m = 2, 
emitting an other photon of the same frequency Ug- If the cavity linewidth n is smaller 
less than 2uJr, only the photons with frequency Ug will survive in the cavity. Since further 
scattering would generate photons with frequencies u — mUr, with m = 3, 5, . . . , which can 
not oscillate in the cavity, then the Hilbert space of the atoms is spanned by only the first 
three recoil momentum levels, m = 0, 1, 2, and Eqs.dH]) and (jlj) reduce to: 

^ = -ga*{t)aci (5) 

= i{6 — ujr)ci + g[a{t)a*co — a*{t)ac2\ (6) 

dc 

= 2i{6 - 2ujr)c2 + ga{t)a*ci (7) 

dt 

— = gNa{t){cocl + C1C2) - na. (8) 

Eqs.(l5])-(IH]) contain fast oscillating terms. They can be eliminated introducing the slowly 
varying variable C2 = C2exp(iAt) and approximating Eqs.Q-® neglecting the fast oscillat- 
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ing terms proportional to exp(±iAt). In this way Eqs.([5])-([8]) reduce to: 

-gaci (9) 



dt 



^ = i{S - ujr)ci + g{a*co - eac2) (10) 
— ^ = 2i{5 — ujr)c2 + gea*Ci (11) 

(JjV 

^ = gN{cQcl + ecicl) - na. (12) 

Eqs.fl9l)- f|T2l) describe the three-level degenerate cascade of the atoms driven by two laser 
fields at frequencies u and uj + 2ujr, respectively, and interacting with the self-generated 
cavity mode at the frequency Ug = uj — uJr- Notice that the second transition rate, from 
m = 1 to m = 2, is proportional to the two pump amplitude ratio, e. 



C. Subradiance in three-level degenerate cascade 

Asymptotically, in a time much longer than I/k, the photons leak the cavity and the 
total polarization in Eq. f|T2|) vanishes: 



Cqc[ + ecic^ = 0. (13) 

On resonance [6 — Ur) and with the atoms initially at rest (co(0) — 1), the variables Cq, Ci, 
C2 and a are real and Eqs.(l9|)- (fT2|) keep invariant the following quantity: 

J = e^cl + cl + 2ecQC2 = e^. (14) 

From it we see that the atoms can not populate completely the final state m = 2 (with 
C2 = 1 and Cq = 0) unless e = 1. Hence, when e ^ 1 the atoms remain in the intermediate 
levels m = and m = 1 in a subradiant state. Condition ( fT3l) . together with the constraints 
f[T^ and the normalization + c\ + c\ = 1, determine univocally the steady-state solution 
reached asymptotically by the atoms. It is easy to show that for e < 1/ "\/3, 



Co — — — 5-, Ci 

1 — 



whereas for e > 1/ v^, 

C0=Y^' Ci = 0, C2 = -^^ (16) 
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Figl2] shows the steady-state populations Pi = |qP for i = 0, 1, 2 plotted vs. e. We observe 
that increasing e from zero the population of the intermediate state, Pi, decreases and the 
population of the final state, P2, increases. They are equal for e = 1/2, with cq = 1/3 
and Ci = — C2 = 2/3. The population of the initial state Pq increases too and reaches a 
local maximum for e = l/VS with cq = 1/2, ci = and C2 = ^3/2. Then, for e > 1/^3 
the intermediate state m = 1 is empty (ci = 0) and the population of the initial state Pq 
decreases to zero for < e < 1; then for e > 1 it increases until it equals the population 

of the final state P2 when e = 1 + ^/2. For e > 1 + -\/2 the population of the initial state, 
Po, is larger than that of the final state, P2. However, this case appears stationary only 
because the semiclassical model neglects spontaneous emission. The quantum treatment, 
reported in the next section, shows that the stationary subradiant state may exist only for 
< e < 1 + 72. 

In order to illustrate how the system evolves toward the subradiance state, figl3] shows 
the time evolution of the field |ap (figlSK), and the three populations, (figlSb) [Pq (dashed 
blue line). Pi (red continuous line) and P2 (dashed-dotted black line)], obtained solving 
the complete equations and (jlj) for gy/N = O.OlcUr, k- = 0.006^;,., 6 = ujr and e = 1/2, 

1.0^ 

W3 

g 

I 0.5- 

Q. 
O 
CL 



0.0 I - . - I — • > ■ , • > ■ 1 

1 2 3 4 5 

8 

FIG. 2: Semiclassical subradiance solution: populations of the three states, Pq (dashed red line), 
Pi (dotted gray line) and P2 ( continuous blue line) as a function of e, as given by Eqs. (fT5]) and 
dUD. 
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with initial condition co(0) = 1, Q^o(O) = and a(0) = 0.01. The final populations are 
Pq = 1/9 and Pi = P2 = 4/9, according to Eq. f|T5l) . In order to test the dependence of the 
subradiance state on the frequency difference A between the two pump fields, figure H] shows 
the asymptotic coherence Ci^2 = l^id^l between the intermediate and the final states vs. A 
for g = O.OlcJr, S = Ur, e = 1/2, k, = O.OOSur (dashed blue line), k, = O-OOGcj^ (dashed-dot red 
line) and k = 0.012co'r (continuous black line). The result shows that subradiance requires a 
very fine tuning of the pump frequency difference near 2uJr, within a precision 5uj <^ gVN- 

As a second example, figs. [5] and [6] show that same case as in figsl3] and |4] but with 
e = 1 + V2. In this case Pi = and Pq = P2 = 1/2. We note that whereas in the case 
e = 1/2 the resonance linewidth of fig S] decreases when the cavity losses k increases, on the 
contrary in the case e = 1 + \/2 the linewidth increases with k and it is about a factor 100 
larger. Hence the subradiance with e = 1/2 is more sensible to the frequency mismatch than 
that with e = 1 + ^/2. 




FIG. 3: Time evolution of the radiated intensity |ap, (a), and the three populations Pq (dashed 
blue line), Pi (red continuous line) and P2 (dashed-dotted black line) vs. t/r, where r = ((7\/iV)~^, 
for gy/N = O.OlcJr, k = OmdiOr, 6 = uor and e = 1/2. 
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III. QUANTUM TREATMENT 



A. The subradiant state 

Let now obtain the subradiant state quantum-mechanically, treating the amphtude c„ as 
bosonic operators c„ with commutation rules [cm,c|^] = 5m,n- Then, according to Eq. f|T3|) 
the subradiant state \sr) satisfies: 



[cqcI + ecicl)\sr) = 



(17) 



(we omit the tilde on 62). It is possible to demonstrate (see Appendix |A]) that for a system 
of atoms (with N even) there are N/2 subradiant states \sr)p, with p = 1,2, .. . , A/2, 
defined as 



\sr 



_ A (-1/26)^ I i2k)\{N - p - ky. 
IP — / . 

fc=0 



k\ \l {p-k)\ 
x\p — k,2k,N — p — k), 



(18) 



where \m,n,l) = \fn) Q\n) 2 and Cp is a normalization constant. The index p is related to 
the population difference between the initial and final states, since Nq — N2 = 2p — A. The 



O 
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FIG. 4: Asymptotic coherence Ci^2 = IC1C2I between the intermediate and the final states vs. 
(A — 2u:r)/g^fN for g = O.Olwr, ^ = w,., e = 1/2, k = COOSu;,. (dashed blue line), k = O.OOGa;^ 
(dashed-dot red line) and k = 0.012LiJr (continuous black line). 
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FIG. 5: Time evolution of the radiated intensity |ap (a) and the three populations Pq (dashed 
blue line), Pi (red continuous line) and P2 (dashed-dotted black line) as a function of t/r, where 
r = (5(\/iV)~^, for the same parameters of figi3]and e = 1 + -v/2- 
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OJ 

0° 0-3- 
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0.1 - 

0.0 I 1 - . . , 

-5 5 

(A-2co^) / (gVN) 




FIG. 6: Asymptotic coherence Co,2 = {cod^l between the initial and the final states vs. (A — 
2uj^)lg\fN for g = O.OluJr, S = u)r, e = 1 + \/2, k = O.OOSuJr (dashed blue line), k = 0.006c<jj. 
(dashed-dot red line) and k = 0.012^;^ (continuous black line). 

case p = corresponds to the state |0,0, A^). The link between the subradiant state \sr)p 
and the semiclassical solution ( fT5l) and ( fT6l) is provided by the correspondence between p 
and the population difference Nq — N2 = N(cl — c^) in the limit ^ 1. For e < l/"\/3, 
p= (A^/2)(l-2e2)/(l-e2) and for e > 1/ VS, p = N[l - Ae^ /{I + e^)^]. As particular cases, 
p = N/3 for e = 1/2 and p = N/A for e = l/-\/3- Furthermore, p = for e = 1 and p = N/2 
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FIG. 7: Quantum subradiance solution: population fraction of the three states, Pq (red circles), 
Pi (black squares) and P2 (blue triangles) vs. e, obtained from Eq. (llSh with N = 32. The dashed 
lines show for comparison the semiclassical solution of figl2j 

for e = 1 + -\/2- Hence 1 + -\/2 is the maximum value of e, giving the following subradiant 
state: 

N/2 



k=0 

A 

X 



A^ 

--k,2k,--k). (19) 



For large A^, the average value of k is {k) = l/4e = {V2 — l)/4, with variance al = {k). 

From the state f[T5]) and the correspondence between p and e we may evaluate the average 
populations Pi = (Ni) /N, with i = 0, 1, 2, as a function of e. The result is compared in fig 13 
with the semiclassical solution (fT5ll and (fT6|) . for A^ = 32. We observe that the quantum 
solution has not a sharp transition at e = l/\/S as the classical one, but there a tail which 
becomes negligible for A^ ^ 1. 
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B. Wigner function 



Here we show the Wigner function of the subradiance state \sr) in order to get some more 
properties of the system. We start from the definition 

H^(«o,«i,«2) = [ e«*"-<^'x(eo,ei,6) , (20) 

where and are complex numbers and x is the characteristic function defined as 

x(eo,ei,6) = {sr\Do{^o)Dii^i)D2msr), (21) 

where Dj{ij) = exp(^jcj — i*jCj) is a displacement operator for the j-th mode. A straight- 
forward calculation, reported in the Appendix [Bl yields 



x(eo,6,e2) = e-(i^°|^+|^^|^+|«^|^)/^x 

p 

Y.Pk V^d^oH L,,{W) L^_,_fc(|6n, (22) 



fc=0 



and 



W(ao,ai,a2) = 0) g-^d^oP+l^iP+I^^P) ^ 
p 

Y,(^k Vfc(4i«or) i:2fc(4i«in L;v-p-fc(4i«2r), (23) 



fc=0 



where 



_ (2fc)!(iV-p-fc)! /ly^ 

and Ln{x) is the Laguerre polynomial. Notice that the Wigner function depends only on 
the modulus of and not from its phase. As expected, in general it is negative due to the 
presence of the Laguerre polynomials. By integrating over the other two mode variables, 
from Eq. fl23|l we obtain the single-mode Wigner functions: 

W{a,) = -e-'\^^\'T{-ir-'/3kL,_,{A\ao\') (25) 

TT ^ — ' 

fe=0 

W{a,) = -e-2|"^l' V/3fcL2fc(4|ai|2) (26) 

71 ^ — ' 

k=0 

2 



W{a,) = le-2|-I^V(-l)^-^-% 

71 ^ — ' 

k=0 

xL^_p_fc(4|a2n. (27) 
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In order to investigate the characteristics of the subradiance state, let's consider some specific 
example. An interesting case is when e = 1/2 and p = N/3, for which the semiclassical theory 
yields Pi = P2 = 4/9. Fig|H](a) shows the probability (3k vs. k for = 36 and p = 12. 
The probability is maximum for k = 8, the average value is (k) = 7.14 and the standard 
deviation is ak = 2.64. The single-mode Wigner functions Wi = W{ai) are shown in figlHl^b- 
d): Wq has a maximum at |ao| =2 and Wi and W2 have a maximum at |ai| = \a2\ =4, 
in agreement with the population values predicted by the semiclassical solution. However, 
Wi differs considerably from W2 with a strong oscillation near = 0, probably due to the 
tail of the distribution at small k observed in figl8](a). The single-mode Wigner functions 
present a pronounced maximum around which they are positive, plus an oscillating quantum 
background. 

As a second example we consider the case e = 1 + ^/2, for which the semiclassical theory 
yields Pi = and Pq = P2 = 1/2. In the quantum model it corresponds to the maximally 
anti-symmetric state f|T9|) with p = N/2. Figj9]J^a) shows the probability 13^ vs. k for A^ = 36 
and p = 18. The probability is maximum for k = and decreases rapidly to zero for larger k, 
with (k) = 0.1 and = 0.35. The single-mode Wigner functions Wi are shown in figin](b-d). 
Wo and W2 are equal and very similar to the Wigner function of the number state \N/2), 



W{a) = (4/7r) exp(— 2|Q;p)LAr/2(4|ap) j2]|. Furthermore, Wi is equal to the vacuum Wigner 
function = (4/7r) exp(— 2|ap). In this case A^/2 pairs of atoms with momentum and 

2hq are produced. 

The three-mode Wigner function fl2^ . after integrating one mode variable, yields the 
following two-mode Wigner functions: 

W{ao,ai) = e-2(l°"l'+l-^l')^(-l)P+'=/5,Vfc(4|aor)L2,(4|ain (28) 

k=0 

W{ao,a2) = e"2(l"°l'+l"^l')f^/5,V,(4|aonL^_p-.(4|a2n (29) 

^ ^ fc=0 

iy(ai,«2) = (^ye"2(l°^l'+l"^l')f^(-l)^"'=/?,V;^(4|aip)L2,(4|a2p) (30) 

fc=0 

Figures [TU| [TT]and[T2] (color online) show the two-mode Wigner functions Wij = W{ai, aj) 
as a function of \ai\ and for i,j = 0, 1, 2 for the case A^ = 36 and p = 12, corresponding 
to e = 1/2 (see figE]): The red color corresponds to a negative value of the functions. We 
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FIG. 8: Subradiant state for = 36, p = 12 and e = 1/2: (a) probability /3k vs. k; (b) Wq vs. 
|q;o|; (c) Wi vs. |ai|; (d) W2 vs. |a2|- 

observe several zones of negativity indicating non classical correlations between the two 
modes. In particular, IVq,! and Wi^2 show strong correlations of the modes and 2 with the 
mode 1, which has strong oscillations, whereas 14^0,2 is larger near the vacuum value (0,0). 
Hence, we can say that qualitatively the modes and 2 behave quasi-classically whereas the 
mode 1 has feature similar to a number state. Figures [T3l shows Wq^2 vs. |ao| and \a2\ for 
the case = 36 and p = 18, corresponding to e = 1 + (see figj9]). The two-mode Wigner 
function looks the product of two single-mode number Wigner functions shown in figjOj^b) 
and (d), with a bi-dimensional regular mesh of positive and negative zones. 



C. Atom statistics 



We calculate now the equal-time intensity correlation and cross-correlation functions, 
defined respectively as: 

(0) = M|M> (31) 

,(^.(0) = AmL, (32) 
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with i = 0,1,2, i j and iVj = c|q. For a classical field there is an upper limit to the 
second-order equal-time cross correlation function given by the Cauchy-Schwartz inequality 



9i:j{0)<[gr'{0)9) 



(33) 



Quantum-mechanical fields, however, can violate this inequality and are instead constrained 
by 



1 



1/2 r 



1 



1 1/2 



(34) 



which reduces to the classical results in the limit of large occupation numbers. We obtain 
the following expressions for the subradiant state: 

^2 



(35) 
(36) 
(37) 



OQ. 






(b) 
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FIG. 9: Subradiant state for N = 36, p = 18 and e = 1 + ^/2■. (a) probability (3^ vs. k; (b) Wq vs. 
|ao|; (c) Wi vs. |ai|; (d) W2 vs. \a2\- 
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FIG. 10: Two-mode Wigner function W{aQ,ai) for N = 36, p = 12. 
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FIG. 11: Two-mode Wigner function W{ao,a2) for N = 36, p = 12. 



where al = {k^) — {ky and 



The cross-correlation functions are 



{k"') = J2k'^Pk 



k=0 



.S(o) = i- 



{No){N^) 



(38) 
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FIG. 12: Two-mode Wigner function W{ai,a2) for N = 36, p = 12. 
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FIG. 13: Two-mode Wigner function W{ao,a2) for N = 3G, p = 18. 
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From ([35])- (HDD it follows that 



(39) 
(40) 
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&0)' = ( g'^\0) + J-) f (0) + T^T V 7^ - 7^ I («) 
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1 







a^iOf = { g?\0) + ^ 1 ( ,f (0) + 7^ 1 - 7^ + 7^ 1 (43) 

showing that glj (0) are consistent with the quantum inequahty ( l34|) . Fig{14] shows the 
intensity correlation functions gl (0) for z = 0,1,2, for the subradiance state (fT8l) . We 
obtain that g2 (0) (continuous red line) is less than unity for all the values of e and g^ (0) 
(dotted blue line) is less than unity for e > 0.56. So anti-bunching occurs for the momentum 
states m = and m = 2, but not for the state m = 1. 

Fig JTSlfTTl show the eventual violation of the Cauchy-Schwartz inequality fl33l) for g^'^j 
(continuous line). The dashed line is the classical upper limit [gl^\0)gf\0)Y/^ and the 
dotted line is the quantum upper limit [(^f^(0) + 1/ {Ni)){gf (0) + l/iNj))]^/^. We found 

(2) (2) (2) 

that gQ { is always consistent with the classical inequality, whereas gQ 2 ^^'i- 9i 2 violate the 
Cauchy-Schwartz inequality (133!) for all the e and for e < 0.47, respectively, showing the 
existence of quantum correlations between the modes m = and m = 2 and between the 
modes m = 1 and m = 2. For e > l/-\/3, (70^2 close to the upper limit of the quantum 
inequality ( !34l) . 



IV. CONCLUSIONS 



We have investigated a possible way to observe subradiance in a Bose-Einstein conden- 
sate in a high-finesse ring cavity, scattering photons from a two-frequency pump laser into 
a single-frequency cavity mode via the quantum collective atomic recoil lasing (QCARL) 
mechanism. Subradiance occurs in a degenerate cascade between three motional levels 
separated by the two-photon momentum recoil hq, where q = k — kg is the momentum 
transfer between pump and cavity mode. The observation of subradiance in momentum 
transitions of cold atomic samples presents several advantages and differences with respect 
to the electronic transitions of excited two-level atoms. First, the momentum transitions 
are not affected by the spontaneous emission if the pump laser is sufficiently detuned from 
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the atomic resonance. Second, the atomic condensates have a long hfe and a very long 
coherence time, allowing the preparation and the further manipulation of the subradiant 
state. Third, the subradiance is realized among collective motional states containing a large 
number of atoms. Subradiance as well superradiance do not need that the dimension of 
the sample is smaller than the radiation wavelength, as in superradiance by excited atoms. 
For these reasons, subradiance between motional states of ultracold atoms may be impor- 
tant for the study of the decoherence-free subspaces sought in quantum information 22|. 



Other recent proposals of realizing subradiance in matter wave require a very fine control 



of single atoms in optical cavities [23|], which can be very problematic experimentally. On 
the other hand, the experimental activity on Superradiant Rayleigh scattering and CARL 
with Bose-Einstein condensates has achieved important progresses and a subradiance 

experiment with BEC in a ring cavity could be feasible with the present day techniques. 
At ultracold temperature and with a coupling constant g^/N much less than the recoil fre- 
quency ujr should be possible, using two laser fields with frequency difference 2uJr, to restrict 
the momentum transition to only the first two recoil momentum states. Recent experiments 



on superradiant scattering from a BEC pumped by a two-frequency laser beam [25|, l26|] have 
shown that the momentum transitions are enhanced by the presence of the second pump 
detuned by 2ujr. However, only inserting the BEC in a high-finesse ring cavity it will pos- 
sible to limit the transition sequence to only two. Then, varying the relative intensity of 
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FIG. 14: Intensity correlation functions vs. e: (0) (dotted blue line), g)^\0) (dashed black line) 
and (0) (continuous red line). 



(2), 
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(2') 

FIG. 15: Cross correlation function gQ[{{S) (continuous red line) vs. e. The dashed black line 
indicates the classical upper limit of the inequality psp : the dotted blue line indicates the quantum 
upper limit of the inequality 
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FIG. 16: Same as in figlTSlbut for g^\{Q). 

the two pump laser beams should be possible to probe the transition from super radiance 
to subradiance. As an example of possible parameters, subradiance could be observed in a 



ring cavity similar to that realized in Tiibingen 27| (with length L = 87mm, beam waist 
w = 100 yum and finesse F = 5 x 10^, about 5 times the presently achieved value) with 
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FIG. 17: Same as in figHSbut for 9^^1(0). 

K = 0.2LUr and g^/N = 0.2LUr. This last value can be obtained using a ^^Rb condensate (with 
Ur = (27r)15kHz) with = 10^ atoms at a temperature of some tens of nK, driven by two 
laser beams with Pq = 174mW, Aq = — (27r)3THz and a frequency difference precision of 
6u! < 5kHz. Subradiance will be reached in about lOO/is, less than the decoherence time of 
the BEG. 



APPENDIX A: THE SUBRADIANCE STATE 

In order to demonstrate Eq. (fT8l) let's consider a state of the form 

N 

\sr) = ^ XI fi^o,n2)\no,ni,n2) 

ni=0 (no,n2) 

where the second sum is over all the pairs (no,?T-2) such that no + ^2 = — rii. When 
substituted in (ITTl) . it yields 
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N-l 



X] f{no,n2)\/no{ni + l)\no - 1, rii + l,n2) 



ni=0 



(no,n2)^no7^0 
"■0+'i2=-'V— "1 



AT 



+ 5Z /('^o,'^2)Vrii(n2 + l)|no,ni - 1,^2 + 1) = 0. (Al) 



ni=l 

{no,n2)^n27^Af 
".o+n.2=-'V— 

After having redefined the indexes ni and n2 in the sums, it becomes 



N 



N-l 



5Z 5Z /('^o + 1,^2)^(^0 + + e X] X] /(no, 7^2 - + 1)?^2 



711=1 



ni=0 



('^0.'^2) 



(™0,n2) 
n2^0 



> |no,ni,n2) = 0. 



(A2) 



The terms in the curl bracket of Eq. (1A2P vanish when the following conditions are met: 

1. In the second sum on rii in the curl bracket the term with rii = is missing: since in 
the first sum no = n2 = when rii = N, it yields /(I, 0) = 0. 

2. In the first sum on Ui in the curl bracket the term with ni = is missing: then the 
second sum on rii with rii = 0, n2 = A^ — no and n2 ^ yields 

/(no,Ar-no-l) = , no = . . . , AT - 1. (A3) 

3. The remaining sum in flA2l) yields: 

^ /(no + 1, n2)\/{no + l)ni + e ^ /(no,n2 - l)\/ini + l)n2 = 0, 



no+n2=N —ni 



(710,712) ^n2^0 
no+n2=N—ni 



(A4) 



with ni = 1, . . . , A^ — 1. In the second sum of Eq. (lA2l) the term n2 = is missing, so 
in the first sum the term with n2 = and uq = N — rii vanishes and yields f{k, 0) = 
with k = 2, . . . , N. Together with the result of item 1, it yields 



f{k,0)=0 , 



AT 



(A5) 
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whereas the remaining terms of Eq. (lA4p yield 

^/{no + l)nif{no + l,N-no-ni) = -e^ (ni + 1)(A^ - no - ni)f{no, iV - % - rii - 1) 

(A6) 

with ni = 1, . . . , — 1 and no = 0, . . . , — ni. 
For ni = 2 Eq.dM]) yields 



v/2(no + l)/(no + l, Ar-no-2) = -e^/3iN - no - 3)/(no, A^-no-3) , no = 0, . . . , A^-S 

(AT) 

Since from Eq.dAlD /(no + 1, A^ - no - 2) = 0, then 

/(no, AT - no - 3) = , no = . . . , A^ - 3. (A8) 

Continuing with all the even values of ni, it is easy to show that 

/(no,A^-no-A;) = , k = 1,3, N - no (odd), no = 0...,N -k. (A9) 

Hence, the only terms different from zero are those with ni = 2g + 1 and q = 0, . . . , N/2 — 1, 
yielding: 



/(no + 1, A^ - no - 2g - 1) = -6^ / ^(^|,^^i)(2° + i) ^^ /(^O' A^ - no - 2g - 2), (AlO) 

where no = 0, . . . , A^ — (2g + 1). Eq. flAlOp provides a recurrence relation for the index no 
with a given n-o + q. In fact the difference between the first and second index of /(a, b) in 
both the left and right terms of Eq.l lAlOl) isA = & — a = A^ — 2(no + q + 1). So, introducing 
the new index p = uq + q + 1, Eq. llAlO l ) can be written, for p = 1, . . . , N/2, as: 



(All) 

By iteration of Eq. (lAlip we obtain 



,f AT ^ I / (2g + l)!! p\iN-p-q)\ 

np-q,N-p-q)=( — ] ,1 (p-,).(Ar-p)/ (P'^-^) ' ^ = 0'---'P 



e 



(A12) 

where (2A;)!! = 2A; ■ (2A; - 2) . . . 2 ■ 1 and (2A; + 1)!! = (2A; + 1) ■ (2A; - 1) . . . 3 ■ 1. Since 
(2g+ l)!!/(2g)!! = (2g)!/(2«g!)2, finally we obtain : 

where Cp = [p\/{N - p)]]^/^ f{p, N - p). 
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APPENDIX B: DERIVATION OF THE SUBRADIANT WIGNER FUNCTION 

We demonstrate Eq. (l22i) . Writing the displacement operator as a product of operators, 
-^(0 = 6xp(~l'CP/2) exp(^c''') exp(— ,^*c) and using the formula 

k 



exp(-rg|»=) = E^\/(^l* 



n=0 



n) 



we obtain 



where 



{k'\D{0\k} = e-'iP'^-Lt(\l;\')Su,, (Bl) 



i.w=i:M)"f*i^ 

n=0 



n 

is the Laguerre Polynomial of order k. Using (IBip in (12T|) with the subradiant state (fTSjl we 
obtain 

p 

j2Pk VfcdeoH i^2fc(i6n ^jv-p-fc(i6n- (B2) 

k=0 

In order to evaluate the Wigner function (!20!) we must calculate an integral of the form: 

U<y) = I d'i e«*°-°*«-l«l'/2L„(|en. (B3) 
Introducing polar coordinates ^ = — exp(i^) and a = |a| exp(2'?/;) it transforms into 

Jo Jo 

POO 

2n / drre-''/'^Lrr,{r^)Jo{2r\a\) (B4) 
Jo 

Where Jo(x) is the Bessel function of zero order. Using the formula 28| : 

dx xe-^'''^^Lnibxy2)Jo{xy) = 



a"+i \2a{b-a) 
we obtain 

27r(-ire-2HY^(4|c,|2) (B6) 
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So, from the definition of Wigner function (120|) and using Eqs. (lB2p and (1B6P we obtain 



1 ^ 

iy(ao,ai,a2) = Ip^k{ao)l2k{ai)lN~p~k{a2) 
n" ^ — ' 



k=0 



which coincides with Eg. (1231) . Using the formula 28 1 

POO 

/ dx e-^/U™(x) = a(l - a)" 
Jo 



we have 



— I d a Im{a) = 1. 



(B7) 



(B8) 



(B9) 



From (EH]), (ET]) and (lB9]) we obtain the expressions i^-l^ and (I28])-(I30]) of the one-mode 
and two-mode reduced Wigner functions, respectively. 
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